Inhomogeneity of the phase space of the damped harmonic oscillator under Levy noise 
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The damped harmonic osciUator under symmetric Levy white noise shows inhomogeneous phase 
space, which is in contrast to the homogeneous one of the same oscillator under the Gaussian white 
noise, as shown in a recent paper [I. M. Sokolov, W. Ebeling, and B. Dybiec, Phys. Rev. E 
83, 041118 (2011)]. The inhomogeneity of the phase space shows certain correlation between the 
coordinate and the velocity of the damped oscillator under symmetric Levy white noise. In the 
present work we further explore the physical origin of these distinguished features and find that 
it is due to the combination of the damped efi'ect and heavy tail of the noise. We demonstrate 
directly this in the reduced coordinate x versus velocity v plots and identify the physics of the 
anti-association of the coordinate and velocity. 

PACS numbers: 05.40.Fb,05.10.Gg,02.50.Ey,82.40.Bj 



I. INTRODUCTION 

The harmonic oscillator is one of several basic models 
that can be exactly solved in textbooks and everyone in 
the physical community is familiar with it. The harmonic 
oscillator can be taken as the most basic model to de- 
scribe the perturbation of an interacting system around 
its equilibrium positions. In reality, the harmonic oscil- 
lators are always under certain environments and their 
influences can be taken into account as the damping of 
the oscillator. However, besides the stable influence of 
the environment, there exists some random perturba- 
tions around the system, which are usually considered 
as the Gaussian white noises, and the dynamics of the 
system is thus physically stochastic. The characteristics 
of the Gaussian white noises are that the noises occur 
in a finite time domain ("white") and the strengths of 
the noises follow a Gaussian distribution. The physics 
of the Gaussian white noises has been intensively stud- 
ied in the past decades. However, in nature many noises 
do not follow the Gaussian distribution but show heavy- 
tailed distributions, for example, a a-stable Levy distri- 
bution [TH5]. Physically, the Levy noise results from the 
strongly collisions between the test particles and/or the 
nonequilibrated heat reservoir abruptly disturbed with a 
local heating to a very high temperature [Gj [7] . Mathe- 
matically, in the Levy noise the distribution of the ran- 
dom variables has a diverging variance. In comparison to 
the conventional Gaussian noises, the heavy-tailed noises 
lead to many novel phenomena observed in some realistic 
systems [MIS. 

In a recent paper [20 , Sokolov et al. investigated 
the phase space property of a damped(including under- 
damped and overdamped cases) harmonic oscillator un- 
der a symmetric Levy white noise. They found that in 
contrast to the homogeneity of the phase space under 
Gaussian white noise, the phase space of the damped 
harmonic oscillator under the general symmetric Levy 
white noises is inhomogeneous, and the inhomogeneous 
extent of the phase space is dependent of a parameter 



a, which describe how far the Levy distribution is away 
from the Gaussian one. The inhomogeneity of the phase 
space shows physically certain dependence between the 
coordinate and the velocity, which is lack in the conven- 
tional Gaussian white noise case. Actually, this prob- 
lem was first addressed and solved in |21j . Their inves- 
tigation showed that the variance of the velocity v{t) of 
the damped oscillator is infinite and imply the particle 
has infinite kinetic energy. As a result, the fluctuation- 
dissipation relation, which can connect the fluctuations 
about equilibrium with the susceptibility of the system to 
external perturbation, is violated in this system. It is the 
breakdown of the fluctuation-dissipation relation in sys- 
tems with heavy-tailed distributions that make studying 
the details of the Levy noise oscillator worth investigating 
in the community. In the present work we revisit this sys- 
tem and further explore the inhomogeneity of the phase 
space. We find that the inhomogeneity can be attributed 
to the interplay between the damped oscillating and the 
Levy noise. To confirm it, we show directly the sample 
points {x, v) obtained numerically in the phase space and 
observe the dependence of the x — v. Here the x and v 
are reduced coordinate and velocity defined by i = x/ax 
and V — v/ay, where and ay are the widths of the dis- 
tributions of the coordinate and velocity [501 121] > which 
are also defined below. The result gives a physically rea- 
sonable interpretation of the inhomogeneity obtained in 
the previous publication pH] . 

The paper is organized as follows. In section [H] we 
present the model and its time-evolution dynamics. In 
section [ni] we discuss the inhomogeneity observed in the 
phase space according time-evolution dynamics. Section 
IV is devoted to a brief summary. 



II. MODEL AND ITS TIME-EVOLUTION 

The problem we are interested in is the phase space 
property of a damped harmonic oscillator under a Levy 
white noise. Here and hereafter, the Levy white noise is 
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FIG. 1. Time-evolution of the Green's functions G{t) and 
Gv{t) for the damped harmonic oscillator with different 
damped cases, (a) The underdamped case {uj = 1,-y — 1), 
(b) the critical case {uj — 1,-y — 2), and (c) the overdamped 
case (a; = 1, 7 = 4). The tfe(fc = 1, 2, ■ • ■ ) in (a) denotes the 
time at which the product of G(t) * Gv{t) changes sign. 



limited to the symmetric discussed in Ref. [20j. 

The model can be described by the Langevin equatfon of 
the coordinate x{t) 

^^{t)+ij^<t) + kx{t)^m. (1) 

where ^(t) denotes the white noise satisfying the symmet- 
ric Levy distribution, k is the elastic constant and 7 is 
the damped coefhcient. For comparison, we use the same 
notations as those in Ref. [20 . It is well-known that 
this model can describe some basic phenomena in the 
nonequilibrium statistical physics. Formally, this equa- 



x{t)^F{t)+ / G{t-t')£,{t')dt', 
Jo 



(2) 
(3) 



where v{t) — ^x{t) denotes the velocity of the oscil- 
lator and F{t) and Fy{t) are decaying functions depen- 
dent of the initial conditions, respectively. Here G{t) and 
[t) are the Green's (response) functions governing the 
time evolution of the system, and they have the relation 
Gv{t) = 4iG{t). It is convenient to take t — t' ~> t', Eqs. 
([2]) and ([3| can be rewritten as 



x{t) - Fit) - 
v{t)=F,{t) 



G{t')at-t')dt', (4) 
Gy{t')C{t-t')dt'. (5) 



By using the Laplace transformation, one can obtain an- 
alytically the expressions of G{t) and Gv{t), which are 
given at different cases. 

Case I: The underdamped case, one has uj — ^/n > 
7/2. Defining fl = ^J>J - i^j^, G{t) and Gy{t) read, 
respectively. 



Git)^'-^^^^^sinm, 
Gy{t) = exp(-7i/2) cos(m) 



7 sin(i7t) 



2VL 



(6) 
(7) 



The typical behaviors of the time evolutions of the G{t) 
and Gy{t) are decaying oscillatingly, as shown in Fig. 

Case II: The critical case, one has uj = ^/k = 7/2, G{t) 
and Gv{t) read, respectively, 

G{t) = texp(-7t/2), (8) 
=exp(-7V2)(l-7V2). (9) 

Case III: The overdamped case, one has uj = ^/k < 
7/2. Defining fl = ^72/4-^2, G{t) and Gy{t) read, 
respectively. 



G(t) = ^^^Pi9^sinh(m), 
G^,(t) = exp(-7i/2) cosh(m) 



7 sinh(ilt) 

m 



(10) 

(11) 



In both the critical and the overdamped cases, G{t) and 
Gv{t) decay without oscillating. These time-evolution 
behaviors of the Green's functions under different cases 
play an important role in understanding the inhomogene- 
ity of the phase space of the damped harmonic oscillator 
under the Levy noise, as discussed later. 
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III. PHASE SPACE PROPERTIES OF THE 
DAMPED OSCILLATOR UNDER LEVY NOISE 

Following Ref . [20] , the characteristic function of the 
joint probability density in the phase space of the damped 
harmonic oscillator at time t is expressed as 



f{k,q,t) = cxp 



\kG{t') + qG^{t')\°'dt' 



(12) 



where a is the Levy index, a is the scale parameter. The 
corresponding joint probability density p{x, v, t) is 



p{x,v,t) 



1 



dk I dqf{k,q,t)exp[—i{kx + qv)] 

) J —OQ 

(13) 

Due to the damping, in long enough time G{t) and 
Gy(t) approach to zero for all three cases, as shown in 
Figjl] and thus the value of f{k,q,t) for given k and q 
reaches a stable value. In this case, it means that the joint 
probability density p{x, v, t) also reaches a fixed point in 
the phase space. 

In the following we solve numerically Eqs. Q and ([5]). 
In the process we consider a Levy noise £,a,D{t) defin- 
ing by jl^^*" £,a,D{t')dt' = La.oi^t), where the value of 
La^oi^t) is given in a Levy process with characteristic 
function as follows 

/CO 
exp{ikLa,D)p{La,D, ^t)dLa^D 
-oo 

exp{ikLc,,D) >= exp(-i:>|fc|"A<).(14) 

Here D = a" denotes the strength of noise. It is 
well-known that La.D has a scale behavior, namely, 
L^ D = -D^/^Xq,! [23!. In Fig. [2] the explicit behav- 
iors of the Levy noise are shown for different indices. For 
the method of generating Levy noise one can refer to Ref. 

[Ml 125]. 

Equations Q and ([s]) are integrated numerically ac- 
cording to the following discrete version 

N-l 



x{NAt) = F{NAt) + {DAt)i ^ G(iAt)CaA^ - <15) 
v{NAt) = Fy{NAt) + {DAt)i ^ G„(iAt)Ca,i(^ -i^) 



where Ca,i(*) is the ith number of the random numbers 
array, which obeys the symmetric Levy distribution with 
index parameter a and scale parameter unity. More de- 
tails on the numerical scheme for integration of stochastic 
differential equations with respect to a-stable noises can 
be found in [U [2311261 [S^ . For simplicity, we just study 
the unity noise, i.e. D = 1. Eqs. (15) and (16) are simu- 
lated by taking initial condition x(0) = 0, viOj = 1, time- 
step At — 0.01 and total N — 2000 steps to guarantee 
time-evolution convergence. We do repeatedly 200000 
times time-evolution of x(t) and v{t) and thus obtain 
200000 sample points in the phase space, which have 
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FIG. 2. Symmetric Levy white noise with different indices a 
(from top to bottom, a = 1.0, 1.5, 1.8, 2.0 ). 



been analyzed statistically. To confirm our numerical 
approach, we reproduce the stationary joint probability 
density p{x,v) for a = 1, as shown in Fig. [3] for the 
underdamped and the overdamped cases. One can note 
that these results are in good agreement with those in 
Ref. [20]. We also obtained the phase distribution p{(j>) 
in Fig. [4|for different Levy indices, which is also consis- 
tent with that presented in Ref. JOj. The angle (p here 
is defined as cf) — arctan('D/5:). 




FIG. 3. The joint probability density p{x,v) for a — 1. The 
results are obtained statistically from 200000 sample points, 
(a) The underdamped case (u = 1,7 = 1) and (b) the over- 
damped case (cj = 1, 7 = 4). 



In Ref. [20] J the inhomogeneity of the joint probabil- 
ity density was attributed to certain association of the 
coordinate and velocity of the harmonic oscillator under 
a Levy noise. Some characteristic features are observed. 
They are that i) the distribution of the coordinate and 
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FIG. 4. Simulated phase distribution p{4i) for different values 
of a. (a) The underdamped case (oi = 1, 7 = 1) and (b) the 
overdamped case (oj = 1 , 7 = 4) . 



velocity is not elliptic, which is in contrast to that un- 
der the Gaussian white noise; ii) the phase angle shows 
a nontrivial distribution; and iii) the dependence of the 
reduced coordinate x — xja^ and the reduced velocity 
V = vja-u is stronger in the overdamped case than in 
the underdamped one. Here Ux = [/q°° |G(i)|"(ii]^/", 
cr„ = |G'«(i)|"fi<]^/". Likewise, the phase angle 
shows more exotic behavior in the overdamped case than 
in the underdamped case. These results were completely 
attributed to the presence of Levy noise. In the following 
we explore carefully the property of the phase space and 
identify the essential physics of these nontrivial phase 
space distributions. 

Without resorting to additional analysis, in Fig. [5] we 
present directly 200000 sample points in the phase space 
with reduced coordinate i, and velocity v of the damped 
oscillator under different white noises. We consider the 
underdamped case (w = 1,7 = 1) in Fig. f5l(a, c) and 
the overdamped case (w = 1,7 4) in Fig. Tsljb, d). For 
the white noises, we consider the Gaussian noise with 
a = 2 and the Levy noise with a = 1.5. Under the Gaus- 
sian white noise, despite of the underdamped case or the 
overdamped one, the distribution of the sample point in 
the phase space is completely homogeneous. This indi- 
cates that the random variables x and v are independent, 
which is consistent with the previous result. However, 
under the Levy noise, the distributions show quite dif- 
ferent features. For the underdamped case, no apparent 
anisotropic distribution is observed but some scattered 
points are far away from {x,v) = (0,0). This result can 
be attributed to the heavy tail of the Levy noise. For the 
overdamped case, the distribution of the sample points 
in the phase space is strongly anisotropic, in particular 
in the regions far away from {x,v) = (0,0). In these re- 
gions, while the sample points can be distributed into the 
regions roughly satisfying x{t) * v{t) > 0, there is no any 
sample points to go into the regions with x{t) * v{t) < 
in the long time limit, which is an interesting observa- 
tion surprised us. This feature is also shown in Fig. 
|6] near {x,v) — (0,0), which is numerically calculated 
by Eqs.(12) and (Il3|). This anisotropic distribution of 
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FIG. 5. Direct distribution of 200000 sample points in the 
phase space for different cases, (a) The underdamped case 
{u — 1,^ — 1) and (b) the overdamped case (a; = 1,7 = 4) 
with Gaussian white noise (a = 2); (c) the underdamped case 
(oj = 1,7 = 1) and (d) the overdamped case (c^; = 1,7 = 4) 
with Levy noise {a = 1.5). The initial condition is taken as 
a;(0) = 0,v{0) — 1, the time-step At = 0.01 and the time- 
evolution steps A'' = 2000. 




FIG. 6. Contour of the joint probability density p{x, v) with 
logrithm scale for Fig. [sj^d) with the same parameters. 



the sample points is just inhomogeneous phase space ob- 
served in Ref. [20]. So, what is the physical origin of this 
inhomogeneity in the phase space? In the following we 
present a detail analysis. 

We check carefully the time-evolution of the Green's 
function G{t) and Gv{t), as shown in Fig. [l] For clar- 
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FIG. 7. (a)-(d) Distribution of the sample points {xk,Vk){k = 
1, • • • , 4) at different time intervals as explained in the context 
and (e) the total distribution of {x, v) after a long enough 
time-evolution for the underdamped oscillator (oj = 1,7 = 1) 
under the Gaussian white noise. The initial condition is taken 
as a;(0) = 0,?;(0) = 1, the time-step At = 0.01 and the time- 
evolution steps — 2000. The time intervals are shown in 
Fig. [Ija). 



ity, we only consider the underdamped case in Fig. [TJ^a) 
with w = 1,7 = 1. When Gv{t) decays with oscillat- 
ing, G{t) increases initially, then also decays with os- 
cillating. So with the time-evolution, the whole time 
domain can be divided into many different intervals in 
which G{t) * Gy{t) changes sign alternatively. For ex- 
ample, when t e {tQ,ti) where to = 0, G{t) > and 
Gy{t) >0, thus G{t)*G^{t) >0. WhentG ^2), G„(t) 
changes sign, Gv{t) < but G{t) is still greater than 
zero, as a result, G{t) * Gy{t) < 0. Furthermore, when 
t € (^2,^3), G{t) * Gy{t) is again greater than zero due to 
that G{t) becomes negative in this time interval. Finally, 
when t e (t3,t4), G{t) * Gy{t) becomes negative since 
G{t) and Gy{t) have different signs in this time interval. 
Of course, the behavior of changing sign of G[t) * G„(t) 




FIG. 8. (a) and (b) Distribution of the sample points 
{xk,Vk){k — 1,2) at different time intervals and (c) the to- 
tal distribution of {x,v) after a long enough time-evolution 
for the overdamped oscillator {cj — 1,7 = 4) under the Gaus- 
sian white noise. The initial condition is taken as a;(0) = 
0,w(0) = 1, the time-step At — 0.01 and the time-evolution 
steps iV = 2000. 



can continue till G{t) and Gy{t) decays away. What does 
this observation mean? Here it is helpful to rewrite Eqs. 
(fT5l) and (fiel) as 



where 



c{NAt)=F{NAt) + {At)^^''Y.Xk, (17) 

k=l 
5 

;(7VAt) =F„(7VAi) + (At)i/"^«fe, (18) 

fe=i 

tk/At 

Xk= G(iAt)Ca,i(iV-*), 

i—tk-i j At 
tfc + i/Ai 

VU= G„(iAt)Ca,l(iV-i), 
i=tfc_i/At 



(19) 
(20) 



where X5 and include all summations for t > t^. In 
Fig. [7] (a)-(d) we show {xi,Vi) for i = 1, • • • ,4 in the 
phase space. It is interesting to find that for {xi,vi) and 
(is, v^) the distribution of the sample points in the phase 
space has a preferred direction along with the first and 
third quadrants, while for {x2,V2) and (£4, {74) the pre- 
ferred direction is dominantly in the second and fourth 
quadrants. This result is easy to understand since in the 
time domains of (ii, ^i) and (is, V3) G{t) and Gv{t) have 
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FIG. 10. (a) and (b) Distribution of the sample points 
{xk,Vk){k — 1,2) at diflerent time intervals and (c) the to- 
tal distribution of (i, v) after a long enough time-evolution 
for the overdamped oscillator (a; = 1,7 = 4) under the Levy 
white noise with a — 1.5. The initial condition is taken as 
x(0) = 0,t;(0) — 1 and the time-step At = 0.01 and the time- 
evolution steps n — 2000. 



FIG. 9. (a)-(d) Distribution of the sample points {xk,Vk){k = 
1, • • • , 4) at different time intervals as explained in the context 
and (e) the total distribution of (x, v) after a long enough 
time-evolution for the underdamped oscillator (tj = 1,7 = 1) 
under the Levy white noise with a — 1.5. The initial condition 
is taken as a;(0) = 0,«(0) = 1, the time-step At = 0.01 and 
the time-evolution steps A'' — 2000. The time intervals are 
shown in Fig. [l|a). 



the same sign and in the time domains of (2:2,^2) and 
{x4,V4) G{t) and Gy{t) have the opposite sign. When all 
sample points are added according to Eqs. ( 17 1 and ( 18 1, 
one obtains Fig. [7][e), which shows a homogeneous dis- 
tribution of the sample points in the phase space of the 
damped oscillator under the Gaussian noise. After un- 
derstanding the homogeneous phase space distribution 
in the underdamped case, we turn into the overdamped 
case with w = 1,7 = 4. The time-evolution of G{t) and 
Gv(t) is presented in Fig. [TJ^c), in which there is no oscil- 
lating. It is noted that only one time moment at which 
Gy{t) changes sign. In this case, strongly anisotropic 
distributions have been observed in Fig. |8][a) and (b). 
However, after merging, one again observes an isotropic 
distribution in the phase space of the damped oscillator 
in the overdamped case under the Gaussian noise. As 



a consequence, despite of the underdamped or the over- 
damped cases, the phase space of the damped oscillator 
is homogeneous if the noise is Gaussian-like. However, if 
the noise is Levy-like, the situation is quite different, as 
shown below. 

In Fig. [9] we consider the phase space distribution 
of the sample points in the underdamped case under 
the Levy white noise. In the different time intervals, 
the preferred directions are similar to those under the 
Gaussian noise, but the difference is that far away from 
{x,v) = (0,0) there are many scattered points to occur, 
which is due to the heavy tail distribution of the Levy 
process. When the sample points in different time in- 
tervals are merged, the distribution of {x,v) is almost 
isotropic, as observed in Fig. |9je). Likewise, some scat- 
tered points are due to the Levy noise. In this case, 
the coordinate and the velocity show weak dependence, 
which is in consistent with the previous work |20j . How- 
ever, the distribution is strongly anisotropic in the over- 
damped shown in Fig. [lOjja) and (b). As a 
consequence, one sees that the distribution of the coordi- 
nate and the velocity of the overdamped oscillator is very 
inhomogeneous under the Levy noise, which is presented 
in Fig. [lO|^c) and also in Fig. [sj^d). In the language 
of Sokolov et al. [20], the coordinate and the velocity 
show strong dependence. According to our analysis, the 
strong (anti) association originates from the interplay be- 
tween the damping and the Levy noise. 
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IV. SUMMARY 



The damped oscillator under the symmetric Levy 
white noise shows inhomogeneous phase space distribu- 
tion as recovered by Sokolov et al. [20, and it was at- 
tributed to the dependence between the coordinate and 
the velocity under the Levy noise. In the present work we 
explore this dependence and find that it can be explained 
by considering the interplay between the damping motion 
of the oscillator and the Levy noise. While the under- 
damped oscillator shows weak anisotropic distribution in 
the phase space, the overdamped oscillator has a strong 



inhomogeneous feature, which is in good agreement with 
the previous result. The essential physics recovered in 
the present work shed a novel light on the dependence 
between the coordinate and the velocity under the Levy 
noise and thus provides a way to understand the physics 
of the damped oscillator under the different noises. 
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